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Abstract. We investigate the second-order von Neumann approach from a 
diagrammatic point-of-view and demonstrate its equivalence with the resonant 
tunneling approximation. Investigation of higher-order diagrams shows that the 
method correctly reproduces the equation of motion for the single-particle reduced 
density matrix of an arbitrary non-interacting many-body system. This explains 
why the method reproduces the current exactly for such systems. Wc go on to 
show, however, that diagrams not included in the method arc needed to calculate 
exactly higher cumulants of the charge transport. This thorough comparison 
sheds light on the validity of all these self consistent second-order approaches. We 
analyze the discrepancy between the noise calculated by our method and the exact 
Levitov formula for a simple non-interacting quantum dot model. Furthermore 
we study the noise of the canyon of current suppression in a two-level dot, a 
phenomenon that requires the inclusion of electron-electron interaction as well as 
higher-order tunneling processes. 
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1. INTRODUCTION 

Understanding transport through interacting nanostructures is of major importance 
for the development of nanoelectronics. The last decades have seen a high activity 
within the field, but calculating the transport through general interacting systems 
remains a challenging task. In contrast, transport through non-interacting systems 
have long-since been well explained by transmission formalism [Tj. In this framework 
the Full Counting Statistics (FCS) was later developed by Lcvitov, Lee, and 
Lesovik [2J. Lately it was also shown how the exact equations of motion (EOM) 
for the reduced density matrix could be derived [31 H] • In a very recent work it was 
shown how this could be used for deriving the FCS in the nonintcracting limit [5]. 
Although very useful in their regime of validity, these methods are restricted to 
non-interacting systems, and can generally not be applied to study the transport of 
confined nanostructures, where electron-electron- interactions are of major importance. 
The transmission formalism can also be used at equilibrium configurations where the 
occupation of the dot is well defined. Here, we are interested in nonequilibrium, where 
the dot is fractionally occupied. 

To deal with the complications of interacting systems, various methods have 
been developed. The most widely used technique is the generalized master equation 
approach which can be derived in many different ways including the real-time 
diagrammatic technique [SI E] by Schoeller, Konig, Schon, and co-workers, as well as 
the Bloch-Redficld approach originally developed in Refs. [HI El LT^ - Comparisons of 
different approaches have been performed in Refs. [TTJ [12] . Bloch-Redficld approaches 
result in a hierarchy of coupled linear differential equations, which must be truncated. 
Keeping only the lowest non- vanishing order in dot-lead couplings one arrives at first- 
order methods, valid in the regime where sequential tunneling is dominant. Then 
one neglects the effects of co-tunneling, giving important contributions for stronger 
couplings [13l HH [15l [16]. Keeping higher-order terms, methods such as the second 
order von Neumann (2vN) approach can be derived |17| . where co-tunneling as well 
as the coherence between the dot states are included. Lately Jin et al. developed a 
method |18j which, at a certain level of approximation was shown to be equivalent to 
the resonant tunneling approximation within the real-time diagrammatic technique. 
Equivalence with the 2vN approach was also suggested. 

The aim of this paper is to investigate the 2vN approach from a diagrammatic 
point-of-view. This allows us to establish a number of exact results for the theory. 
Firstly, we demonstrate the equivalence between the 2vN method and the resonant 
tunneling approximation, (and thus the method of Jin et al). This proof is enabled by 
the compact notation of Liouville diagrams. We then consider non-interacting systems 
and show that, although the 2vN approach yields the correct EOM for the reduced 
single-particle density matrix p, enabling an exact calculation of the current |18j , the 
EOM for the complete reduced many-body density matrix p is not exact. This is shown 
via the cancelling of higher order diagrams. To see how the cancelling works we have 
apply Keldysh diagrams. These provide a more detailed description, and the Liouville 
diagrams are actually sums over distinct Keldysh diagrams. Building on the work of 
Ref. |19j we consider the 2vN approach with the inclusion of counting fields. Using 
our diagrammatic method, we are able to explain why the shot-noise calculated for a 
non-interacting system is not exact, even though the current is. If fact, we find that 
the inclusion of diagrams of order T 4 is required to reproduce the noise correctly (r is 



Diagrams for equations of motion, current, and noise for the 2vN approach 



3 



the lead-dot coupling energy, which contains the square of the tunneling elements ). 
Finally, we consider the noise calculated numerically within the 2vN approach for two 
illustrative models. We compare our noise results for the single-resonant level (SRL) 
model with the exact results and demonstrate that, despite the limitations discussed 
above, the 2vN approach still gives a good description of the noise properties of this 
model for T < ksT. We then investigate the noise and Fano-factor for the two-level 
spinlcss system where a canyon of current suppression was previously observed |20| 12 1 j . 

The remainder of the paper is organized as follows: In Sec. [2] and Appendix A we 
present the 2vN approach with the inclusion of counting fields Sec. [3] presents 

the proof, using Liouvillian perturbation theory, of the equivalence between the 2vN 
approach and the resonant tunneling approximation within the real-time diagrammatic 
technique. Sec. 0] gives a brief overview of the real-time diagrams employed in this 
work. Using diagrammatic techniques, we first show in Sec. [S] that the 2vN approach 
gives the correct EOM for the reduced density matrix of a single-level system coupled 
to reservoirs, but that this result does not generalize to the reduced density matrix p 
of a multi-level many-body system. However, the EOM for the reduced single-particle 
density matrix is shown to exactly reproduced for non-interacting many-body systems, 
which allows for exact calculation of the current. In Sec. [6] it is shown that this does 
not hold for the noise, which cannot be reproduced exactly even in the non-interacting 
limit. Sections [7] and [5] contain the results of our 2vN calculations for the single- level 
system and the interacting two- level system. We then conclude in section [9] 

2. The second order von Neumann method with counting fields 

A common transport Hamiltonian describes transitions (T) between two leads via a 
quantum dot (D) and can be written as 

H = HLoads+HD+HT- (1) 

We assume that the leads are Fermi liquids of electrons with dispersion Em- With I 
we distinguish the different leads. Index k is the (quasi-)momentum quantum number, 
but it could also contain additional properties like spin: 

%Leads = "^^EktC^Ckl- (2) 

ki 

We describe the main system, i.e. the quantum dot, in its diagonal basis using the 
states | a) which can be many-particle states including also bosonic degrees of freedom: 

n D = ^2E a \a){a\. (3) 

a 

Responsible for transport are tunneling processes between the main system and the 
leads. These processes are described by terms where the state of the main system 
changes while at the same time an electron leaves or enters a reservoir. The tunneling 
Hamiltonian including counting fields, Xe with I signifying at which lead electrons are 
counted, is given by [35] 

«t(A) = Tba(ki)\b)(a\c ke e~ lXe ' 2 + h.c. (4) 

k£,ab 

X When we refer to the order of a method, we refer to the power of T, and not that of the tunneling 
elements. 
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In this section we use the convention that dot state \b) contains one more electron than 
dot state a), and so on. In this formalism one can derive the generalized Liouville- 
von-Neumann equation 

i jLp( A) t) = H + (X)p(X, t) - p(X, t)U~ (A), (5) 

where, like in the rest of the paper, we set K = 1 for the sake of simplicity. We have 
introduced the A-dependcnt Hamiltonian 

W^A) = ^Lcad S +H D +^T(±A). (6) 

In Laplace-space Eq. ((2) reads 

izp(X, z) = ip(X, 0) + H + (X)p(X, z) - p(X, z)U-(X). (7) 

For states of the whole system we choose a basis of tensor products \ag) = |a)® \g), 
where |a) describes the complete many particle state of the main system and \g) 
describes the many particle state of the leads. 

We introduce the following quantities 

Wb'b = y^/Vq;fe fl , 
9 

4>ba{kl) = ^p hg _ U . ag , (8) 
9 

Xab(k£) = ^2 Pag;bg-ki 
9 

where g — k£ denotes a lead state where an electron with quantum number k in lead £ 
has been removed. Note that Xab(k£) ^ <t>t a (k£) due to the inclusion of counting hclds. 
Wb'b are the elements of the reduced density matrix. In the EOM these elements couple 
to cj> and x which are the coherences arising from the superpositions of an electron in 
the leads and the dot. These in turn couple to elements where two electrons are moved 
between leads and dot. To break this hierarchy of equations the EOM is truncated as 
described below. This results in a closed EOM which can be solved numerically. 
In order to achieve this we use the following approximations: 

• Factorization and Thermodynamic limit: Wherever an occupation operator 
of a lead state appears we replace it by a Fermi function fki of the respective lead. 
For example: 

^2{ag\cuCke4,, t ,p(X,z)\bg) 

9 

« fu^2(ag\cl, e ,p(X,z)\bg). (9) 

9 

This approximation is valid under the assumption that relaxation in the leads 
is quick, compared to the tunneling of individual fc-states, so that equilibrium is 
restored between each tunneling event. This is well justified in the limit of an 
infinte number of lead states. Hence the occupations on the leads will not change 
due to the tunneling of single electrons, even not for large times. 
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• Truncation: We neglect all matrix elements where the bra- and the ket-state 
differ by more than two single-particle states of the leads. For example: 

(ag\c kiei cl 2£2 c k3 i 3 p(X,z)\bg) « 0, 

if kiii, k-iii and £^£3 are different. (10) 

This represents neglection of higher order tunneling events, and is not valid in 
the regime T ^s> ksT. 

Following the derivations in Refs. [T71 Qj5] we can write down the EOM for <fi, x, 
and w, see Appendix A. Previously, implementations of this scheme have solved 
these equations numerically [TTJ [T5] (and see section [7]) but here we are interested 
in determining analytical properties of the approach. For this purpose the real-time 
diagrammatic technique [6] [7] is more suitable, as discussed in Sections [4] and [5] 



3. The 2vN method in Liouville space 

The structure of the 2vN scheme is most easily seen using the diagrams of 
Liouvillian perturbation theory (LPT) [23] [24] [25] . In particular, using this approach 
we demonstrate the equivalence of the 2vN scheme with the resonant-tunneling 
approximation. This approximation was introduced for the SRL in Ref. |26j and 
is defined in diagrammatic terms as retaining all irreducible diagrams in the kernel 
where a vertical cut crosses at most two lead contractions. Here we derive an explicit 
expression for the self-energy of the 2vN approach, which shows that the 2vN approach 
corresponds to exactly this approximation for an arbitrary system. We first discuss 
the theory without counting field, and include it at the end. 



3.1. Liouville space 

We begin by giving the essentials of LPT, but the reader is referred to Refs. [2"3I |2~4"1 |2"T1 
|2"S] for full details. The von Neumann equation for the evolution of the total density 
matrix under Hamiltonian Eq. (JTJ) reads: 

p(t) = -i[H,p(t)}=r,p(t). (11) 

which defines the Liouvillian super-operator C. In accordance with the decomposition 
of the Hamiltonian, C consists of three parts: C = C rcs + Cd + Ct with £ los = 
— i [%Lcads,» ]) £d = —i[Hn, 9 ], and Ct = — i [Ht, • ]• To ease book-keeping, we 
introduce a compact single index "1" to denote the triplet of indices (£i,fci,£i). The 
notation T refers to the triple(— £1, k\, l\). The first index £1 = ± describes whether a 
reservoir operator is a creation or annihilation operator: 

ai = a (lklil = ! C L' ^ 1=+ ■ (12) 
I CMn ?i = - 

Similarly we define 

91 I E ab T ba (kl)\b)(a\, a = - ' 1 J 

such that the tunnel Hamiltonian can be written Ht = ^\a\g\, with all summations 
left implicit. The sign £ enters here because the reservoir operator always comes first 
in this expression for Ht, unlike in Eq. ([4]). 
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The tunnel Liouvillian £t = ~~ i [Ci a i3i; •] consists of system and bath operators 
acting from both the left and the right, i.e. on both Keldysh branches. We introduce 
the superscript Keldysh index p = ± to describe these two possibilities and define 
corresponding super-operators in Liouvillc space. For the reservoir, we define £-space 
super-operator A via its action on the density operator p: 

At\p))»{ aiP > P = + . (14) 
[ pax, p = - 

Following Rcf. [21], we define the system super-operators G via 

G p O = a? x ( P = + . (15) 

t P = ~ 

To avoid confusion we explicitly state the G is thus not a Green's function. The object 
a p is a dot-space super-operator with matrix elements 

I v\ x x J 1 ' N s -N s ,= even 

KW = W^| p) N s -N s ,= odd ' (16) 

where, N s is the number of electrons in state s. The tunnel Liouvillian can then be 
written 

Ct = -itxpaVAlG*. (17) 
3.2. The 2vN EOM in Liouvillc space 

Arranging the elements Wb'b into the vector \w)) and elements <pba{kih) into vector 
|0(1))) the 2vN EOM (without counting field), Eq. (fA~T|) and Eq. (fO|) . read in 
Liouvillc space 

z\w)) - \w))(t Q ) = C D \w)) + £t(1)|#1)» (i 8 ) 

i 

= (*i+£d) +2>(1)H 

+ 5^,(1,4)^(4), (19) 

4 

where, T(l) = iGf . We have also defined T/(l) = £f> 2 /(-&P2Wi)Gf with Fermi 
function /(w) = [e' 3 " + l] -1 , and finally the block 

M>(1,4) = -P3/(-€3P3W 3 )x (20) 
5l4Gf " ^Gf - <5l3^42Gf — Gf 

° z — X3 — x\ — Lb * z — X4 — x\ — Lb 

where summation over index 2 and 3 is implied, and where X\ = — + Pi) 

with wi the energy of lead mode {k\,l\) and p\ the chemical potential of lead l\. 
Diagrammatically, M/(l,4) is represented by: 



G 3 G 3 G 2 G 3 



where symbol G represents a tunnel vertex, the lines between the G:s represent 
free propagation, and the lines on top indicate contraction of lead operators. 
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3.3. The 2vN self- energy 

An expression for the system self-energy within the 2vN method can be derived by 
iterating Eq. (fH)]) . Defining the free system propagator 0(1) = [z — x\ — £d] _1 we 
obtain 

|0(1))} = 0(1) |r y (l)H + £ M/(l, 2)|0(2))) | 

= 0(1)T/(1)H + 0(1)5^M/(1, 2) jfi(2)2>(2)|«;»+ 0(2) ^M f (2, 3)|0(3)))1 



= n(l)T / (l)H + 0(1)M/(1, 2)0(2)T/(2)H 
+ 0(1)M/(1, 2)0(2)M/(2, 3)0(3)T/(3)H + . . 

Substituting into Eq. (fT8)) and solving gives 



(21) 



K*)» = r 1 w > (*o))) (22) 

z — £d — S(z) 

with the (non-Markovian) self-energy 

£(z) = T(1)0(1)7>(1) + T(1)0(1)M/(1, 2)0(2)T/(2) 

+ T(1)0(1)M/(1, 2)0(2)M/(2, 3)0(3)7/(3) + . . . (23) 

with summation over all indices implied. This we can write in a compact form. Let 
T be a vector with elements T(l), T/ be a vector with elements 1/(1). Furthermore 
let M be the matrix with elements M/(l, 2) and finally, O be a diagonal matrix with 
elements fi(l)fli 2- The self-energy can then be written 

£(z) = T- n-T/ + T- n- M- n-T/ 

+ TOMOMOT/ + ..., (24) 
which can clearly be resummed to give 

E(z) = T — — — — • O • Tf (25) 

w 1 — n ■ M 1 v ' 

which is a very nice clean result. 

This self-energy can be expanded in terms of the diagrams of Liouvillc 

perturbation theory, see Refs. f23j [28j [25]. For example, the product of MO • M 

entering in Eq. (|23|) formally corresponds to the diagrammatic expression, 



I I X 
G — G- \ G — G G — G 



+ 



X 



G — G- 



l I + I — 

G — G G — G 



where all indices have been suppressed. When diagrams are multiplied unpaired 
contraction lines are connected. As each diagram in M only has one unpaired 
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contraction line it is ensured that there will never be more than two contraction 
lines at any time. By including the terminations T and TV in Eq. (|23[) and remaining 
free propagators, one arrives at the diagrammatic expression for the self energy 

i 1 I I i F=i I 

£(z) =G-G+G-G-G-G+G-G-G-G 

I FF FF I I Pi 1 

~\~ G — G — G — G — G — G ~\~ G — G — G — G — G — G 

I ff i I ff — i 

+ G-G-G-G-G-G + G-G-G-G-G-G + ..., (26) 

From this expansion it is clear that the 2vN self-energy contains the set of all diagrams 
in which there are at most two contraction lines passing over any given point. This 
exactly proves the equivalence between the 2vN approach and the resonant tunneling 
approximation. 

Finally we define the propagator in Laplace space 

p{u) = n(w)p(t ) = . 1 , . p(t ), (27) 

— iuj + Lb + Zj{—ilj) 

which will be used in Sec. 6. 



3.4- Counting statistics 

Counting fields can easily be added in the kernel Eq. (|25p by replacing the bare tunnel 
vertices with their gauge-transformed counterparts: G Pl — > G\ x e l ^ lPlXe ^^ 2 [28]. In 
the exponent here, £i corresponds to whether an electron is created or annihilated 
(see just before Eq. (fT2")l ) and p\ is the Keldysh index. With this replacement, the 
non-Markovian A-dependent 2vN self-energy then reads 

S (A;,) = T(A)- i _ n 1 M(A) .n-T / (A) (28) 

with matrices 

T(A;1) = iG p { e i?lPlAi i /2 

T>(A;1) = ip2f(-ZiP2Ui)G? e^P^'i/ 2 

M(A;1,4)= (29) 

- P^f(-^P^) {<^i4GfO(l + 3)Gf e l ^-P^/ 2 

-5 13 6±2Gfn(l + ^)Gfe- i ^ p ^ ,2 e i ^ x ^' 2 \ 

where e.g. 11(1 + 4) denotes the free propagator [z — x\ — — £d] _1 - 

Derivatives of S(A; z) with respect to A then give the FCS. The current and 
finite-frequency noise can be written in a "current block" notation as HH] 

(I) = e«J lQ (0))), (30) 
e 2 

S((jj) = —{{J2a{u)p s t + Jla(u))Il((j)J lb (cj))) 

- 2tt5(uj)(I) 2 + terms with (w -> -ui), (31) 

where ((...)) = denotes expectation with respect to the stationary state of the system 
and where J\ a , Jif,, and J^b are current blocks (see Appendix |Appendix B| for explicit 
forms). The diagrams contributing to Ji a (u) and Ji&(w) are equivalent to those 
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contributing to the self energy E(z = —ioS) except that one tunneling vertex has 
been replaced by a current vertex. In J\ a (uj) the leftmost vertex has been replaced, 
while in Ji b (u>) any replacement is possible. Such a replacement changes the value of 
a diagram by a factor of 1/2 and furthermore introduces a sign change if the current 
vertex is moved to the other branch. We also define J2a(w) where two tunneling 
vertices have been replaced by current vertices. One of the current vertices must be 
to the left, the other can be at any position. The reason that J\ a and J2a needs a 
current vertex at the leftmost position is that when the trace is taken and the Keldysh 
contour closed all diagrams with a tunneling vertex to the left will cancel each other. 

4. Diagrammatic description 

Our aim is now to investigate when and why the 2vN approach is exact. For this 
purpose the LPT used in the previous section is too compact. We can, however, 
use the result of the previous section that the 2vN approximation is equivalent to 
the resonant tunneling approximation within the real time diagrammatic technique. 
From here on we will therefore only use two-branch Keldysh diagrams (which the LPT 
diagrams are summations over). The upper and lower contour of the Keldysh diagram 
give the time evolution of the bra- and ket-states respectively This corresponds to 
different ordering on the contours so that the contour ordering of the tunneling events 
agrees with the direction of time along the upper branch, while it is opposite to the 
direction of time on the lower branch. One can view this as the time running backwards 
along the lower branch. This section is devoted to a discussion on how the EOM for 
the reduced density matrix of a central quantum system coupled to external leads can 
be derived using such diagrammatic methods. 

In deriving the EOM we follow the diagrammatic notation of Ref. (T2], where 
more details on how the diagrams are evaluated can be found. Assuming that the 
total density matrix is a direct product of the initial states of the dot and leads at the 
time to, when the couplings between leads and dot is switched on, the EOM for the 
elements of the reduced density matrix p b -y can be written as a closed set of linear 
equations 

dt = -i{E b - E b < )p b . v 

+ ]T f drK^'(t-r)p a . a! (r), (32) 

aa' J *° 

where the first term on the right corresponds to the unitary evolution of the dot 
disregarding couplings to leads and the second term is generated by tunneling events. 
Unlike in Sec. [5] it is not assumed that state \b) contains one more electron than state 
| a). Here the coefficients K^, can be calculated using Keldysh diagrams as outlined 
in Ref. [12]. All possible irreducible diagrams connecting the states \b), \b') on the 
left with the states \a), \a') on the right contribute to the coefficient Kfig, . In an 
irreducible diagram any vertical cut crosses at least one lead contraction arrow. We 
use the following notation \b) = \a, /3, ...), where anti-symmetrization of the many- 
particle state is implicitly assumed, to specify which single-particle states |a), 
that compose the many-body state \b). This way of writing the many-particle state 
as product states is generally not possible for interacting systems but can always be 
done in the non-interacting case that we focus on. 
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Here we give a brief summary of the diagrammatic rules, using an example 
diagram giving the time evolution on /9i,2,3 ; i,2,3 m terms of see Fig. [TJ More 
detailed discussions can be found in Refs. [SJEIIT^]. Throughout the paper we use the 
convention that time increases from right to left. 



1,2,3 1,2 1 




Figure 1. Example of a 3:rd order diagram included in the 2vN method. 

1) Each vertex where a lead contraction arrow ends has a related tunneling 
probability Tf, a (kl), representing that the dot state changes from many-particle state 
| a) to many-particle state \b) by the annihilation of an electron in lead-state k in lead 
I. Arrows pointing out of a branch correspond to Hermitean conjugated processes. 

2) In a contraction, the lead operator with the larger time argument comes first 
unless the right vertex of the contraction is on the lower branch. For Fig.[T]this results 
in the Fermi factors: f{Ek) [1 — f(Ek>)] f(Ek><), where we have included the lead index 
I in k. 

3) The factors in the denominator of the diagram are given by the dashed vertical 
lines of Fig. [T] Each line gives a factor which is given by (the energy on the lower 
contour) - (the energy on the upper contour) + (energies of arrows pointing right) - 
(energies of arrows pointing left). Furthermore i0 + should be added to each factor. 

4) The sign of the diagram is given by (— 1)-"+™=+"* x (Fermion sign), where N is 
the order of the diagram in T, n c is the number of crossings of contraction lines, and 
ni is the number of vertices on the lower contour. In Fig.[T]^V = 3, n c = 2, and n\ = 4 
so that the first sign factor becomes —1. Comparing with Ref. [12] the Fermion sign 
results from the ordering of the dot operators. For the diagram in Fig. Q] the Fermion 
sign is given by 010103020302 11) = — 11), i.e. the Fermion sign is negative, resulting 
in an overall positive sign in front of the diagram. 

5. Cancellation of diagrams for non-interacting systems 

In this Section we show how all diagrams discarded in the resonant tunneling 
approximation, i.e. the diagrams not included in the 2vN approach, cancel each other 
in the EOM for the single-particle reduced density matrix in the non-interacting limit. 
The single-particle reduced density matrix is defined as 

P K u = Tr dot [a^p] , (33) 

where the trace is taken over the dot system, p and p are related to each other in 
the following way: The diagonal states of p give the probability of the corresponding 
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many-body state being occupied. This probability corresponds to that exactly the 
single-particle states that compose the many body state are occupied, the others are 
empty. The diagonal terms of p give the probability that the corresponding single- 
particle state is occupied, without considering the occupations of the other single- 
particle states. 

The exact EOM for the single-particle reduced density matrix of non-interacting 
systems is a major strength of the 2vN approach as important physical observablcs 
such as the current can be exactly calculated from this knowledge. 

For single level QDs a canceling partner diagram can always be found for each 
diagram by changing the branch of the second and third leftmost vertices in the 
diagram, see Figs. [2] and [3] For multi-level QDs the diagrams can either be canceled 
by methods similar to those used for the single level QDs, or be grouped into families, 
see Fig. [5] where all diagrams have the same absolute value but half contribute to the 
EOM with plus signs and half with minus signs. 

5.1. Single resonant level 

To explain on a diagrammatic basis why the resonant tunneling approximation gives 
the right EOM for the single-particle reduced density matrix and the right current 
we consider diagrams with three dot-lead excitations at the same time. We begin by 
noting that diagrams where different vertices are contracted necessarily have different 
denominators. When trying to find out how the diagrams cancel a good strategy is 
therefore to look for diagrams where the same vertices are contracted, i.e. diagrams 
that correspond to the same Liouvillian diagram (see Sec. [3]). Since it is a single- level 
quantum dot, we are limited to occupations of |0) or |1). This means that at every 
second vertex we must create an electron in the dot while the other vertices must 
annihilate an electron. Specifically this means that if we create at vertex 4 we must 
do the same at vertex 3 if they are on different branches, and if we create at 3 we 
must annihilate at 4 if they are on the same branch. 

A canceling diagram can be found for any third-order diagram not included in the 
2vN approach by changing branch of vertex 3 and 4. This results in the same change 
of occupation at the two branches so that the energy denominators are unchanged. 
As we only moved the left vertices in the contraction, the contributions from the 
lead contractions (Fermi factors) are also unchanged. Since we are dealing with a 
single-particle system, we can neglect the Fcrmion sign and the total sign is given 
by (— 1)-™+™':+"'!. Clearly N is unchanged. The change in ni is or ±2, i.e. an even 
number. The change in n c is given by the number of dot operators that are commuted. 
The two moved vertices will both be commuted with the unmoved left vertex 5. Then 
there is the additional commutation between the two moved vertices. In total n c is 
therefore an odd number. 

We will illustrate the canceling using an example diagram seen in Fig. [2] The 
red numbers on gray backgrounds give the state of the quantum dot, while the black 
numbers label the vertices. Changing branch of vertex 3 and 4 for the diagram in 
Fig. [2] results in the cancelling partner diagram shown in Fig. [3] 

For higher-order diagrams the reasoning is much the same. A canceling diagram 
can again be found by moving the second and third leftmost vertices to the other 
branch. In the same way as for the third-order diagrams the change in number 
of crossings is an odd number so that the diagrams cancel. This explains on a 
diagrammatic basis why the 2vN method gives the exact EOM for p and the exact 



Diagrams for equations of motion, current, and noise for the 2vN approach 12 

1 4 2 1 
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Figure 2. Example of a 3:rd order diagram, for a single level system, that is not 
included in the 2vN approach. The red numbers on gray backgrounds denote the 
dot occupation while the black numbers label the indices. 



1 3 2 1 




Figure 3. Canceling partner of the diagram in Fig[2] 

current for the single level |17j . 

The above described method explains how to cancel diagrams where three dot- 
lead excitations exist at the same time, i.e. the type of diagrams not included in the 
2vN approach. We emphasize that higher order diagrams where the maximum number 
of dot-lead excitations are restricted to two, see e.g. Fig. 21 cannot be canceled in this 
way. Changing branch of the same vertices as before, in Fig. |4] labeled by 2 and 4 
since the vertices have been shifted in time, changes the number of crossings in the 
diagram with an even number. As a result the two diagrams do not cancel. This 
demonstrates the improvement of the 2vN approach compared to pure second order 
approaches where diagrams such as Fig. H] are not included. 

5.2. Canceling of diagrams for multi-level systems 

For multiple level QDs it is generally not possible to cancel the diagrams in the 
above described way. The reason is that there is no requirement that electrons 
are created/annihilated at every second vertex as many different states in the dot 
can be active. In order for the diagrams to cancel the two moved vertices must 
create/annihilate the same state in the dot. Generally it is not possible to find 
two such vertices in many-body diagrams. As a result, the 2vN approach does not 
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Figure 4. Example of a higher-order diagram included in the 2vN approach that 
is not canceled by the above described method. 

give the correct EOM for p for many-body systems even in the noninteracting limit. 
Importantly, we will see that the EOM for the single-particle reduced density matrix 
is still correct in this limit. 

Following the definition of the single-particle reduced density matrix p we have 
for non- interacting systems, where the many-particle states always can be written as 
direct products of single-particle states: 

Pp;v = ^2 Pl*+c;u+c- (34) 
c: (n,v(£c) 

Here, the state /Lt+c is the many-particle state c with a particle added in single-particle 
state p. To get the time evolution of p^ v in terms of p a - a ' we must therefore sum over 
all diagrams giving the time evolution of p^+ c -.v+c in terms of p a:a '. In this sum many 
of the diagrams will have the same energy denominators and Fermi factors, i.e. the 
values of the diagrams are equal, which means that they cancel if they have opposite 
sign. 

Before continuing we introduce the concept of free vertices: The free vertices of 
a diagrams are those among the N left vertices where states other than p and v are 
active. To explain how the canceling works we divide all diagrams of third order and 
higher into four groups: 

1) No state occurs more than once among the free vertices. 

2) Some states occur twice among the free vertices. 

3) Some states occur three times or more among the free vertices. 

4) There are no free vertices. This can only happen if // ^ v. 

It is clear that any diagram falls into one of theses groups. Below we explain how 
the canceling works for each group. 

Group 1: Wc will see that a family of diagrams can be constructed by changing 
branch of the left vertices where states other than p or v are active. It is clear that 
every diagram belongs to one and only one family. Moving the vertices in such a way 
does not affect the energy denominators of the diagram, and since only the left vertices 
are moved the contributions from the Fermi-factors are also unaffected. The diagrams 
within a family will thus cancel if equal numbers contribute with plus and minus sign. 
It should be noted that the energy denominators of the different diagrams in a family 
are only equal in the non-interacting limit. For interacting systems the diagrams do 
not cancel. 
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The number of free vertices will be denoted by n and the number of these that 
are on the upper branch will be denoted by k. The number of diagrams in a family 
with k out of the n free vertices on the upper branch is then given by (^) . 




1,3 1,2,3 1,2 1 1,2 1,2,3 1,2 1 



d) i o 




1 1,3 1,2,3 1,2 1 



Figure 5. A family of canceling diagrams in group 1. 

We first show that all diagrams with a given k have the same sign. In the 
expression (— l^+^+ii the number of crossings between contraction lines changes 
when we move one vertex from the upper to the lower and one from the lower to the 
upper. However, this is canceled by the additional change in the Fermion sign, as 
each change in the number of crossings is related to the commutation of two Fermion 
operators of different dot states. In total there will therefore not be a sign change. 
Next we show that the sign changes when k is changed by one. Moving a vertex can 
result in a change in the number of crossing contraction lines, but this number is again 
the same as the number of fermion operators that is commuted due to the change. 
Left is only the change in sign due to that the number of vertices on the lower branch 
has changed by one. In the general case one therefore gets a sum 

E("l) fc (T) =0 forn>l. (35) 

k=0 ^ ' 

Since each diagram belongs to one and only one family and each family gives no 
contribution we conclude that this group of diagrams gives no contribution. 

An example family of diagrams belonging to group 1 is shown in Fig. [5] Here, 
\fi) = \v) = |1), and on the right we have the state |0). This family of diagrams thus 
contribute to the time evolution of px;i m terms of po;0- F° r this family of diagrams 
n = 2, corresponding to the states |2) and 3). 

Group 2: When some states occurs twice among the free vertices only the leftmost 
of these two vertices can be moved to the other branch. If the right of these vertices 
is moved, a state will be created or annihilated twice in a row. However, all diagrams 
in this group have at least one vertex that can be moved to the other branch which 
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results in a family of canceling diagrams according to Eq. (|35[) . An example of such a 
diagram is shown in Fig. [S] Moving the right vertex marked with a red square results 
in that state \2) is annihilated twice in a row. 




1,2 

Figure 6. A family of canceling diagrams in group 2. Only the vertex marked 
with a green square can be moved to the other branch. 

Group 3: If the same state |a) occurs three times or more among the free vertices 
the diagram can be canceled by exactly the method used for single-level diagrams. 
When moving the vertices corresponding to state |a) crossings with contraction lines 
for other states will occur. However, this sign change is canceled by the commutation 
of Fermion operators as described above. Vertices and contraction lines not belonging 
to \a) can therefore be neglected which results in a single level diagram in state a. 

Group 4: Diagrams with no free vertices can for N > 3 only occur if n ^ v. If 
there are no free vertices in the diagram it means that no states other than pi or v 
occur in the diagram. Thus at least two of the N vertices to the left must act on 
the same state. These two vertices must necessarily be on the same branch and one 
creates the state while the other annihilates the state. If there is no creation vertex 
of the same state on the other branch to the left of these two vertices, a canceling 
diagram is found by changing branch of the two vertices, see Fig. [7] a). In the other 
case one changes branch on the two creation vertices, see Fig. UJb). Fig. [7] only shows 
the two possibilities for the three left vertices as the cancelling of the diagrams can be 
explained without looking at the right part of the diagram. Moving the vertices like 
this works in the same way as for single level diagrams. It results in two diagrams 
with opposite sign but equal Fermi-factors and energy denominators. 

To conclude one can divide all diagrams into those that cancel directly for p, and 
those that cancel only for p. The ones that cancel directly for p can be canceled by 
moving two vertices in a way that does not change b and b'. To cancel the rest of the 
diagrams for p the diagrams are divided into families. All diagrams within a family 
contributes to the time evolution of the same element in p. All diagrams in the family 
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Figure 7. Demonstration of how diagrams belonging to group 4, i.e. diagrams 
with no free vertices, are canceled. 



have the same value but half contribute with a plus sign and half with a minus sign. 
As a result of the exact EOM for p, the current can be calculated exactly for any 
non- interacting system, as outlined in Ref. [18]. 

It is interesting to note that although the 2vN method docs not correctly 
reproduce p, it can in some cases still be calculated exactly using the knowledge that 
the EOM for p is exact. For systems where the probability to tunnel into superpositions 
of the dot states can be neglected, such as the non-interacting Anderson model or 
dots where the level spacing greatly exceeds the coupling strengths so that the secular 
approximation can be performed, each dot level can be treated independently from 
the others. This allows us to write down the following expression for the diagonal 
elements of p, 

Pb:b = Y[Pa\a ~PP;P)> ( 36 ) 

aeb P0 

where the first product is over all single-particle states composing b and the second 
product is over those not belonging to b. For such systems the 2vN method thus allows 
for an exact derivation of the EOM also for p. We illustrate how this works for the 
non-interacting Anderson model. Here the elements of p can be calculated from p as 

P0;0 = (1 -P-htX 1 -Vi;\), 

Ptit =Pt;t( 1 -Pl;l)> 

Pl-A = C 1 -Pt;t)PUh ( 37 ) 
Pd-,d = Pt,fPi;h 

where p±d denotes the doubly occupied state. We again remark that this is only 
possible for non-interacting systems where the occupations of the different single 
particle levels are independent. 
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6. Calculating the noise for single resonant level systems 

To investigate if the noise is exact we will study the constituting parts of Eq. (|3T|) . 
From the above discussion it is evident that the time evolution operator in Laplace 
space II(w), Eq. (|27p . is correctly reproduced by the 2vN approach for single-particle 
systems. Furthermore the diagrams in the J\ a block of Eq. (f3Tj) excluded from the 
2vN scheme can be canceled by the same method as for single-particle systems. 



4 2 




5 3 1 

Figure 8. Example of a 3:rd order noise diagram for J2 a which cannot be canceled 
by the above described method. The green dots correspond to current vertices 
while the black represent normal tunneling vertices. The trace in Eq. J 3 1 1 ) is 
represented by the closing of the Keldysh contour. 



3 2 10 




7 6 5 4 

Figure 9. Example of a 4:th order noise diagram for J2 a which cannot be canceled 
by the above described method. The green dots correspond to current vertices 
while the black represent normal tunneling vertices. The trace in Eq. H31H is 
represented by the closing of the Keldysh contour. 

However, the current vertex in blocks Jib and Ji a that occurs at positions other 
than leftmost prevents the canceling of diagrams. Moving such a vertex to the other 
branch results in an additional sign change so that the diagrams add upp instead of 
cancelling. This is illustrated in Fig. [8] which cannot be canceled by moving vertex 3. 
Thus, this 3:rd order diagram is likely to contribute to the noise. 

There are also 4:th order diagrams that contribute to the shot noise, see e.g. 
Fig. [S 
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It should be noted that although the counting statistics was introduced for finite 
frequency noise, the results also holds at zero frequency so that the noise cannot be 
correctly calculated by the 2vN approach in this limit. 

Going to fifth order or higher, one can always find two normal tunnel vertices 
that can be moved without resulting in a diagram with double occupation. 4:th order 
in r is thus required to calculate the shot noise exactly for the single level. Each order 
of cumulants in the FCS adds an additional current vertex. As a result one needs to 
include two more orders of diagrams, which quickly results in great numerical efforts 
being required. 

It should be pointed out that we have not presented a proof that the 3:rd and 4:th 
order noise diagrams do not cancel. Rather we have shown that the method by which 
the diagrams for the EOM cancel, does not work for the noise diagrams. It could be 
argued that the noise diagrams cancel by some other manner, but previous analytical 
work have shown that the 2vN method does not reproduce the noise exactly for the 
single resonant level [T5]. That 4:th order diagrams are required is also expected 
from the transmission formula, see Ref. pQ. In the low-temperature limit the noise 
is proportional to T(l — T), T being the transmission function. Here the T 2 term is 
proportional to T 4 . It is interesting to note that although the 2vN approach does not 
give the correct noise, it can still be obtained by inserting the correct transmission 
function from the 2vN method into the transmission formula. In the non-interacting 
limit the transmission function can be obtained from the current as it does not depend 
on the applied voltage. 

7. The single resonant level 

Having seen that higher-order terms are required to correctly reproduce the noise of 
single-particle systems it is of interest to see how large the discrepancy is between the 
2vN method and the exact results using the transmission formula. 

To get an impression of this discrepancy, we apply the 2vN approximation to the 
noise in the single resonant level model. 

Assuming the (constant) tunneling rates T a = 2ir J^ fe \tk a \ S(E — E^ a ) and 
T = T L +Tr, and defining the quantity B e (E) = 2irJ2k *fc£0io(z, A; M)5{E — Eu), we 
obtain the EOM for the single resonant level 




izu>n(z, A) 



iwn(t 





(38) 
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These equations have simple analytical solutions at infinite bias, when the tunneling 
rates are constant or have a Lorentzian shape [T5]. At finite bias however, we need a 
numerical solution and hence must transform the function B(E) into a mapping on 
a discrete and finite set. Two numerical solution procedures (A and B) have been 
independently developed. 

• A Discrete energies: In our first approach [17] we choose a set of equidistant 
energies and approximate the integrals with sums over these energies. This 
requires a cutoff at some energy and is therefore most efficient for scenarios where 
a sufficiently small energy interval determines the transport properties. 

• B Cutoff in residues: One can show that in the crucial complex half plane the 
quantity B(E) has the same poles as the Fermi function [30]. By this knowledge 
the integrals can be approximated with the residues that are closest to the real 
axis. This leads to an inaccuracy in the value of the Fermi functions. Thus, this 
method is most efficient when the transport happens at energies where the left 
and the right Fermi function differ sufficiently. 

The two methods also differed in the way the noise was evaluated: 

• A Numerical time evolution: We evaluate the time evolution of the cumulant 
generating function S(X,t) = — In [%2 b Wbb(^,t)]. After a short time, S(\,t) 
becomes linear in time and its slope s(A) can be determined. The current and 
noise can then be calculated from ^ and §v§ respectively. 

• B Long time expansion in Laplace space: By comparing the factors in 
front of the system states w a b in Eq. ()38[) we can extract a non-Markovian master 
equation. As soon as this is known, we can use a Taylor expansion around 2 = 
and the techniques described in Refs. [3TJ [35] to calculate the noise. 

In approach A the EOM is solved in time space by performing the Markov 
approximation for the density matrix elements containing two dot-lead excitations. 
In Laplace space this corresponds to replacing all z to the right of the equality 
sign in Eq. (|A.2[) with positive infinitesimals while keeping the other z in the EOM. 
Using somewhat hand-waving arguments it is easy to realize that doing the Markov 
approximation at a higher-order is less of an approximation. The z to the right of 
the equality sign in Eq. (|A.2[) contribute when the rest of the denominator is small. 
This region can be defined as a volume in the iV-dimensional energy space, N being 
the order at which the Markov approximation is done. As N increases this volume 
becomes smaller compared to the entire integration volume, which reduces the effect 
of neglecting the z. 

The accuracy of the 2vN approach and the effect of the partial Markovian 
approximation is investigated in Fig. [10] where the noise, i.e. the second cumulant 
((J 2 )), is shown as a function of the bias, Vbias, for two different coupling strengths. 
The dot level is positioned in the middle of the bias window and we assume equal 
couplings to left and right lead, = Tr. Partial Markovian results from approach 
A are compared with the non-Markovian results of approach B, and the exact non- 
interacting results of the transmission formula. For both coupling strengths it can be 
seen that the effect of performing the partial Markov approximation as described 
above, is very small. Partial Markovian calculations were also performed using 
approach B. This gave results indistinguishable from approach A. For weak couplings 
such as Fig. [TU] a) the 2vN results agree very well with the exact transmission 
results as higher-order terms are of less importance (note the y-axis scale). For the 
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stronger couplings of Fig. [TU] b) the agreement is good in the low bias limit. As the 
bias is increased the phase space for higher-order processes increases which causes a 
discrepancy between the 2vN results and the exact transmission formula. 
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Figure 10. The noise {(l 2 )) for a single-level model for a) Fl = Fr = O.lfcsT, 
b) = Tfi = kgT, and the dot level positioned at the middle of the bias window. 
In the strong coupling case the ncglcction of higher order tunneling results in 
a discrepancy between the 2vN results and the exact transmission formula as 
expected. 



8. Noise and Fano factor for the canyon of current suppression 

So far the success and shortcoming of the 2vN method have been studied. In this 
section we change focus and use the partial Markovian version A of the method to 
calculate the noise and Fano factor of the canyon of current suppression, previously 
investigated experimentally [5D| and theoretically [3TJ [33J GUI [33 [Ml [37] . The canyon 
appears as a suppression of current in both sequential and co-tunneling regimes, close 
to degeneracy in two-level spinless systems. 

In the single-particle eigenbasis of the dot, the system Hamiltonian is given by: 

H = H dot + i?icads + H T , (39) 

H dot =E 1 d\d 1 +E 2 4d 2 + Ud\di4d2, (40) 

-ffleads = ^ EkC^Ckl, (41) 

k,e=L/n. 

Ht = (tiid\+t t2 dl)c M + R.c, (42) 

k,e=L/R 

where we have assumed that the couplings t^ are independent of fc and Tu{E) = 
2Trtf i p with a constant density of states po(E) = J^k ^(^fc — E), for —D < Ef, < D. 
In the simulations a large bandwidth D is used, assuming wide conduction bands of 
the leads. The operators di (d\) and cm (cL) are annihilation (creation) operators of 
electrons in the dot and leads, respectively. In Eq. PH]) the charging energy U is due 
to Coulomb repulsion between the electrons when both dot states are filled. 
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Noise calculations for this system have previously been performed in Rcf. |38j 
under the assumption that one level was very weakly coupled to the reservoirs. Here 
we report results for the two level system without this assumption, demonstrating 
the versatility of the 2vN approach. Noise calculations have also been performed for 
similar systems such as serial double quantum dots [39], albeit in the weak coupling 
regime. 

We parametrize the energy levels as 
AE 

E 1/2 = ±— -E g - U/2, (43) 

where AE is the splitting between the two levels and E g is a common shift of the 
levels. 

We study couplings of the type 

tn = t, t m = t, t L2 = -at, t m = at, (44) 

where the asymmetry parameter a is chosen to be real, as such couplings were shown 
to be essential for observation of the canyon [21] . 

Figures llla )-c) show normalized current ((I)) /Vbias, normalized noise 
((/ 2 )) /Vbia S , and the Fano factor F = {{i 2 }) /((/)) as a function of AE and E g 
using the parameters Vbias = 15k bT, T L i = T m = k B T, T L2 = T R2 = k B T/4, and 
U = 25k B T corresponding to the parameters of Fig. 2 (e) in Ref. [21] . 

Fig. [TT] a) shows the canyon of current suppression discussed in Ref. [21] . To 
understand the next order cumulant of charge transport we look at the results for 
the Fano factor. The most prominent feature in the plot are the two hills located at 
AE/k B T = -10, E g /k B T = 5 and AE/k B T = 10, E g /k B T = -5. These hills result 
form the trapping of the electron in a low conducting state. The level configuration 
at AE/k B T = 10, E g /k B T = —5, marked by a white cross in Fig. [TT]c), is shown in 
Fig. [IT] d). Here, state |2) of the dot is at most times filled, which corresponds to a 
very low conductance, as current is carried by co-tunneling processes through cither 
the weakly coupled state |2) or through state |1) which is far from the bias window due 
to Coulomb interaction. However, sometimes inelastic co-tunneling events empty |2) 
and fill |1), resulting in a highly conductive system. As a result the electron transport 
is bunched which corresponds to as large noise signal and a high Fano factor. The 
simultaneous treatment of sequential tunneling and co-tunneling in the presence of 
Coulomb interaction is therefore essential for describing the charge transport in this 
regime. Thus, the 2vN approach is ideal for this type of calculations. Neglecting 
the effects of finite temperature and level broadening this phenomenon occurs at level 
configurations 

Vbias , „ V h 

< hj\ < — - — , 



2 

Vbias 

2 



< E x + U, 



E 2 < (45) 

and similar for the hill located at AE/k B T = —10, E g /k B T = 5. For the parameters 
considered in Fig.[TT]the limits are set by — ¥&bb. < E\ and E 2 < — ^ fl§ -, corresponding 
to the areas restricted by the red lines in Fig. [TT]c). Furthermore, the level splitting 
must be smaller than the bias to enable conservation of energy. This requirement does 
not impose any restrictions in Fig.QT]c) as < Vbias in this plot. For \ AE\ > Vbi as 
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AE/k B T 

Figure 11. a) The normalized current ((/)) /Vbias, b) the normalized noise 
((J 2 )) /Vbias, in units of e 2 /h, and c) the Fano factor F = ((I 2 )) / ((/». The 
parameters are given by V bias = 15k B T, T L1 = T R1 = k B T, T L2 = T fl2 = 
fcgT/4, and U = 25fcsT. d) The level positions at the white cross in Fig. c). 

the current drastically decreases, and the Fano factor drops to a value close to unity 
as the inelastic co-tunneling processes no longer are possible, see Fig. [T2] 

This phenomenon is different from the dynamical channel blockade described in 
Refs. [501I3T], where both levels are located within the bias window and tunneling can 
be described within the sequential regime. The regime considered in the current work 
gives a stronger bunching of electrons as the dot at most times is in a state where there 
is no state inside the bias window that can carry the electron. Indeed, in the regime 
where both states are in the bias window corresponding to the dynamical channel 
blockade, a Fano factor close to one is observed for the studied parameters. There is 
an exception to this is close to degeneracy, AE = 0. Here the Fano factor exhibits a 
"ridge" resulting from interference between the two dot states. When the dot states 
are degenerate, electrons can tunnel into any superposition of these states. Especially 
it can happen that an electron tunnels into a linear combination with zero coupling to 
the right lead. When this happens the electron is trapped in the dot [32] as it cannot 
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Figure 12. The Fano factor as a function of AE for E g = —5ksT. The other 
parameters are given by Vbias = 15fc s T. F L1 = F R1 = k B T, F L2 = V R2 = 
kgT/A, and U = 25kgT. For AE > Vbias inelastic co-tunneling processes are no 
longer possible, resulting in a lower Fano factor. 

easily tunnel back to the left lead due to the high bias. This explain the canyon of 
current suppression in the regime where both levels are well inside the bias window. 
Furthermore it results in strong bunching of charge transport [43], i.e. the ridge in 
the Fano factor, as rare co-tunncling events allow the trapped electron to escape back 
to the left lead. Unlike dynamical channel blockade, this is not a modulation of the 
charge transport where the weakly coupled state controls the transport through the 
strongly coupled state. Instead it is an interference effect between the two levels. 

For ~5ksT < E g < 5ksT and AE = 0, there is no state in the bias window. 
Even if the dot leaves the non-conducting state, its conductance is still low as current 
can only be carried by co-tunneling processes. Thus, the ridge is weaker in this regime. 

Estimating the error in these simulations is somewhat more difficult as it is not 
possible to benchmark against the exact transmission formula for interacting systems. 
However, the results presented in Ref. [30] show that at couplings of T ~ kBT and 
Vbias = 15fcsT the noise agrees very well with the transmission formula for the single- 
level dot. This suggests that we are in a regime where higher order tunneling events 
are of less importance and the 2vN method can be trusted to provide accurate results. 

9. Conclusions 

The 2vN method and the resonant tunneling approximation within the real-time 
diagrammatic approach were compared using Liouvillian perturbation theory. This 
established the equivalence between the two methods and the recently developed 
method of Ref. [TS]. Using diagrammatic techniques it was shown how diagrams 
of third and higher order in the lead-dot coupling energy, not included in the 2vN 
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approach, canceled in the EOM for the single-particle reduced density matrix of non- 
interacting systems. This enables the exact calculation of the current. However, it was 
seen that to correctly reproduce the EOM for the many-body density matrix, or the 
Full Counting Statistics, diagrams of higher-order were needed. The discrepancy in the 
noise between the 2vN method and the exact transmission formula was investigated 
in the non-interacting limit. For stronger lead-dot couplings an increased discrepancy 
was observed at larger bias, due to an increased phase space for higher-order tunneling. 
Finally, the noise of the canyon of current suppression was calculated using the 2vN 
method. While the current and noise showed a canyon around level degeneracy, the 
Fano factor exhibited a ridge, as well as local maxima due to electron bunching. 
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Appendix A. Equation of motion for the 2vN approach including 
counting fields 

Using the notation of Sec. [2] the EOM, originally derived in Ref. [19], is given by 
(z - E b + E b ,)w w = +Y, eiXe/2 (Tba{ke) X ab>{ke) - (j> ba {kt)T^ a {M)) 

+ J2 e" lAf/2 {T* b (kl)4, cW (kt) - X bc(ke)T cb ,(k£)) 

c,kt 
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(z - E c + E b + E k )(/) cb (k£) = 

+ Y, e iAf/2 T cb ,(k£)f ke w b , b - fktWcc>T c , b (k£) e- iX 'l 2 

V c' 

[e iX * T cb ,(k'e')f kU '(l> b ' a (k£) - f* kl i>c V {k'£')T b , a {k£)e^~^l 2 \ T* a (k'£') 



E 



E k + E k > - (E c - E a ) + z 

a,b',k'i' y ' 

[r k ,r<j>cb'(k£)T b , a (k'£i) - e KV+^)/2 T cb ,(ke)f M cj )b , a (k'e')] T* a {k't) 
+ £t„ E k +E k , -(E c -E a ) + z 

a,b,k'£ f 

v T cb ,(k'£>) [f k 'i4b'a{k£)T; a {k'e)e^' -e^'+ x ^' 2 T b , a {k£)f klX a b {k' £')] 
+ trt, e , E k - E k , - {E b , -E b ) + z 

v T cV (k'£') [T* lbf (k'l')f^ t ,<l> c , b (k£) - r kl Xb'Ak'£')T c , b {k£)e^~^/ 2 ] 
v }%„ E k -E k ,-{E b ,-E b) + z 

v [fk'i-<t>cb'{k£)T* b ,{k'£') - e '(^-V)/2 T cb ,{k£)f MX b>c{k'£')\ T db (k< £') 
vfc* E k -E k> -(E c -E c ,) + z 

J e -i V T* dc {k'£')f* k , t 4 dc ,{k£) - f^ X cd(k'£')T dc ,(k£) e -i(^+Ai)/aj Tc , b (k'£') 

+ E k -E k ,-{E c -E d ) + z 

T* dc {k'£') [r k ,^ d Ak£)T c , b {k'£') e~ iX * - e'^-V)^ T dd (k£)f kl ^ b (k' £')\ 
+ £t„ E k + E k , - (E d -E b ) + z 

v T* dc {k'£') [T dc ,{k'£')f k , t< p c , b {k£) - r M <p dc >{k'£')T c , b {k£)e-^'+^/ 2 ] 
£rt„ E k + E k , -(E d ~E b ) + z 

(A.2) 

here = 1 — f k e. The EOM for \ can be obtained by taking the complex conjugate 
of Eq. IA.2I followed by replacing Xe with — Xi, 4>l a with Xab> and x* b with 4> ba . It 
should be noted that 4>, Xi an d w au depend on the counting fields A^. 

Appendix B. LPT current blocks 

The blocks used in the current and noise expressions, Eqs. (ptU|) and (f5T|) . can be 
derived from the following two primitive blocks: 



jW(z 1 ,z ) = G-X' + X'-G (B.l) 



I i=i I I i=i 1 

X'-G-G-G + G-X'-G-G 

Z Zq Z 2i Z Z 

I r=^ 1 I F=i 1 

G-G-X'-G + G-G-G-X' 

Z\ Z\ Zq Z\ Z\ Zi 
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i jW(z 2 , Zl , z ) = X' - X' + ±X" -G+^G-X" 

Z\ Zq Z2 



-X' -X' -G-G 

Zi Zq Zq 



-X' - G-G - X' 

21 21 21 



1 

X 


i — 

' -G- 


1 

X' 


1 

- G 




21 21 






1 

G 


1 

-X' - 


1 

X' 


1 

- G 




22 21 




zo 



I I 1 I I I 1 I 

-G — X' — G — X' +G-G-X' -X' (B.2) 

22 Zl Zi Z2 Z2 Zi 



I i=i — I I i=i — I 

V -G-G-G +±G-X"-G-G 

ZQ Zq Zq ' Z2 Zq Zq 



r 



— G — X" — G ~^\G — G — G — X" -(-..., 

2 2 22 2 ' 2 2 2 2 2 2 

where we have written out explicitly only the sequential and direct co-tunneling 
contributions. In the diagrams, X' and X" represent the first and second derivatives 
of the A-dependent tunnel vertex evaluated at A = 0. The blocks of Eq. (f5(J)) and 
Eq. (pTTj) are obtained as 

Jio(w) =jW(-iu,-iu); 

J 2a (w) = ^J' (2) (-*w J -««; I 0) (B.3) 

where in the latter two blocks, we can throw away all diagrams with a leftmost G- 
or X"-vcrtex due to the occurrence of these blocks to the far left in the current and 
noise expressions [2"5] , 
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